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Abstract 

The longitudinal asymmetry due to Z° exchange is calculated in quasi-elastic 
electron-deuteron scattering at momentum transfers \Q 2 \ ~ 0.1 GeV 2 relevant for 
the SAMPLE experiment. The deuteron and pn scattering-state wave functions are 
obtained from solutions of a Schrodinger equation with the Argonne v±s potential. 
Electromagnetic and weak neutral one- and two-nucleon currents are included in 
the calculation. The two-nucleon currents of pion range are shown to be identical 
to those derived in Chiral Perturbation Theory. The results indicate that two-body 
contributions to the asymmetry are small (~ 0.2%) around the quasi-elastic peak, 
but become relatively more significant (~ 3%) in the high-energy wing of the quasi- 
elastic peak. 



pacs: 13.10.+q, 21.45,+v, 25.30.Fj 



1 Introduction 



With results from the SAMPLE experiment available 0, there is now great interest in 
parity-violating electron-deuteron scattering. The SAMPLE Collaboration measures at 
Bates the longitudinal asymmetry in polarized quasi-elastic electron scattering on the 
proton [|IJ and deuteron 0. The asymmetries are sensitive to the nucleon's strange mag- 
netic and axial form factors. If other effects are under control, the two measurements 
can be used to determine the values of the form factors at \Q 2 \ = 0.1 GeV 2 . Neglecting 
two-nucleon current effects, preliminary results are in disagreement with theoretical es- 
timates, in particular for the axial contribution ||. The purpose of the work presented 
here is to compute a theoretical value for the deuteron asymmetry including two-nucleon 
currents. 

Neutral charge and current one-nucleon operators are well known (see, for example, 
Ref. ||). The contribution of such operators to the asymmetry was computed at various 
momentum transfers in Ref. || . Several theoretical issues were studied in detail. For the 
kinematical region relevant to SAMPLE, final-state interactions were found to be impor- 
tant. It was also found that (except for very low momentum transfers) the asymmetry 
in the vicinity of the quasi-elastic peak is fairly independent of the choice of two-nucleon 
(NN) potential. 

Two-nucleon charge and current operators have also been studied to some extent. 
Heavy-meson exchange contributions were considered in Ref. ||. They were shown to 
be unimportant in a calculation of the asymmetry that neglected final-state interactions. 
In Ref. @, an impulse approximation modified to incorporate gauge invariance was 
employed. The effects of parity-violating NN interactions on the deuteron wave function 
were found to be small. Pion-exchange currents were included in the computation of the 
asymmetry, but only in electromagnetic sector. 

In the present work we extend these earlier calculations. The leading currents are 
calculated in Chiral Perturbation Theory (ChPT) || [| supplemented by a successful 
phenomenological model ||10|| . The one-nucleon currents considered here include phe- 
nomenological form factors and have the same form as those in Ref. || . We also consider 
pion-exchange currents in both electromagnetic and weak contributions, as well as effects 



from heavier mesons, evaluated using the Riska prescription [11|. The asymmetry is cal- 



culated with deuteron and final-state wave functions obtained from a realistic potential, 



the Argonne v±s model [12 



We present results for the kinematical region of interest to SAMPLE. (Our calculation 
can be repeated at other momentum transfers, such as those of the planned SAMPLE- 



Lite experiment |L3[ and JLab's GO experiment ||14j| .) We verify that when two-nucleon 



effects are turned off, our results are in agreement with those of Ref. ||. We study the 
effects of two-body currents on the asymmetry, both near the quasi-elastic peak where 
one-body processes should dominate, and away from the quasi-elastic peak where these 
two-body currents could be important. We are able to confirm that most of the two- 
nucleon contributions to the asymmetry are due to currents of pion range. Our results 
show that near the quasi-elastic peak two-body currents give a small contribution to 
the asymmetry. Away from the peak, they become more important, and can increase 
the magnitude of the asymmetry by as much as 3%. We find that the contribution to 
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the asymmetry associated with the electromagnetic-axial current interference response 
function is about 20%. The overall effect of two-nucleon currents on the data of the 
SAMPLE experiment is indeed small but not negligible, and is being incorporated in the 
data analysis ||. 

A summary of the relevant formulas for the calculation of the asymmetry is given in 
Sec. H Section |3.1| contains the one- and two-body electromagnetic currents, and Sec. 



contains the one- and two-body neutral weak currents. The connection between our 
phenomenological approach and ChPT can be found in Sec. f|, while the numerical com- 
putation of the asymmetry is described in Sec. |5|. Results and conclusions are presented 
in Sec. |]. 



2 The Parity- Violating Asymmetry 

Parity-violating electron-nucleus scattering results from interference of amplitudes asso- 
ciated with photon and Z° exchanges, shown in Fig. |I|. The initial and final electron 
(nucleus) four-momenta are labelled by and k 1 ^ (P^ and P'* 1 ), respectively, while the 
four-momentum transfer is defined as Q M = k^ — k' 11 = (u, q). The amplitudes for the 
processes in Fig. [I] are then given by (T5| 

M = --^(ikT + M z ) , (1) 
M-> = uYujl fi , (2) 

= A7T ^2 ~ u 7 \9v + 9A ls) u 3vji , I 3 ) 

where a and are the fine-structure constant and Fermi constant for muon decay, 
respectively, gy 1 = — 1 + 4sin 2 6 ) n/ and = 1 are the Standard Model values for the 
neutral-current couplings to the electron given in terms of the Weinberg angle 8w, u and 
u' are the initial and final electron spinors, and jjf and j^' CT denote matrix elements of 
the electromagnetic and weak neutral currents, i.e. 

jjf = (f\f'°(0M = (P}M) , ft(q)) , (4) 

and similarly for j'tf- Here \i) and |/) are the initial and final nuclear states. Note that 
in the amplitude M z the Q 2 dependence of the Z° propagator has been ignored, since 
here we restricted ourselves to \Q 2 \ m|. 

The parity- violating asymmetry in the quasi-elastic regime is given by 

/ da + da~ \ I ( da + do~ \ , . 

" \d£ldu ~ dttdu J I [dnduj + dVtdu J ' 



where da h /dQdu is the inclusive cross section for scattering of an incident electron with 
helicity h — ±1. It is easily seen that, to leading order, 

A oc \ ^ ,„ ; . 6 
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Figure 1: (a) Single photon- and (b) Z°-exchange diagrams in parity- violating quasi-elastic 
electron-deuteron scattering. The blobs denote the nuclear currents. 



Standard manipulations then lead to the following expression for the asymmetry in the 



extreme relativistic limit for the electron [jnj 

12 „(e)„,_ p7,0 , „(eV,_ d7,0 , Je), p7,5 



.4 



2^2 a vlRY 1 + v T RY f 

where the v's are defined in terms of electron kinematical variables, 



(7) 



Q 4 



v T = tan 2 (0/ '2) + 



2q 2 



tan(#/2Wtan 2 (#/2) 



\Q 2 



(9) 
(10) 



9 being the electron scattering angle in the laboratory, while the i?'s are the nuclear 
electro-weak response functions, which depend on q and u, to be defined below. To 
this end, it is first convenient to separate the weak current j z ' u into its vector j°' a and 
axial-vector j 5,cr components, and to write correspondingly 



■Z,a 

hi 



■0,<T 1 -5,(7 

Jfi +3fi 



(p°/ i (q)J?.(q)) + (p 5 /i (q)J5 i (q)) • 
The response functions can then be expressed as 



i f 



R^ a (q } w ) = £ ]T 5( W + m, - E f) ® 

i f 



[ID 



:i2) 



(13) 
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lT/(q, u:) = Y,T,^ + m i- E f) 3 

i f 

where m{ is the mass of the target (assumed at rest in the laboratory), Ef is the energy 
of the final nuclear state (in general, a scattering state), and in Eqs. ( |T2| ) and ( |TB| ) the 
superscript a is either 7 or 0. Note that there is a sum over the final states and an average 
over the initial spin projection states of the target, as implied by the notation m the 
expressions above for the i?'s, it has been assumed that the three-momentum transfer q 
is along the z-axis, which defines the spin quantization axis for the nuclear states. 

The model for the nuclear electro-weak currents is discussed in the next sections, while 
the calculation of the deuteron response functions is described in Sec. §. 



(14) 



3 Electro- Weak Charge and Current Operators 
3.1 Electromagnetic Operators 

The nuclear charge and current operators consist of one- and two-body terms that operate 
on the nucleon degrees of freedom: 



p 7 (q) = EpT 1 (q) + E4 ! (q) > 

i i<j 

j 7 (q) = EjJ' 1 (q)+Ej5 2 (q) • 



(15) 
(16) 



The one-body operators p/' 1 and j/' 1 have the standard expressions obtained from a rela- 
tivistic reduction of the covariant single- nucleon current, and are listed below for conve- 
nience. The charge operator is written as 



p 7,1 (q) = p 7 ,NR(q) + p 7 ,R C (q) , 



(17) 



with 



p 7 ic(q) 



1 



1 & e iqr * - 



4m 2 



x Pi )e iq ^ , (19) 



'l + \Q 2 \/Am 2 

where \Q 2 \ = q 2 — uo 2 > is the four- momentum transfer defined earlier, and m is the 
nucleon mass. The current operator is expressed as 

j 7ll (q) = [p, , e^] + -^qx cr L , (20) 

where [• ■ • , ■ • ■]+ denotes the anticommutator. The following definitions have been intro- 
duced: 



1 r 



e,; = - 



GU\Q 2 \) + G V E (\Q 2 \)T; 



2 

1 r 



[Li 



GU\Q Z \) + G v M (\Q 2 \)r z , 



(21) 
(22) 
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and p, <r, and r are the nucleon's momentum, Pauli spin and isospin operators, respec- 
tively. The two terms proportional to 1/ m 2 in p]'^ c are the well known Darwin- Foldy and 
spin-orbit relativistic corrections [jnj, respectively. The dipole parametrization is used 
for the isoscalar (S) and isovector (V) combinations of the electric and magnetic nucleon 
form factors (including the Galster form for the electric neutron form factor 



The most important features of the two-body parts of the electromagnetic current 
operator are summarized below. The reader is referred to Refs. [[T^, [H], [HJ for a derivation 
and listing of their explicit expressions. 



3.1.1 Two-body current operators 

The two-body current operator has "model-independent" and "model-dependent" com- 
ponents, in the classification scheme of Riska |TT| . The model-independent terms are 



obtained from the two-nucleon interaction (in the present study the Argonne Vi$ inter- 
action is employed), and by construction satisfy current conservation with it. The 
leading operator is the isovector "7r-like" current obtained from the isospin-dependent 
spin-spin and tensor interactions. The latter also generate an isovector "p-like" current, 
while additional model-independent isoscalar and isovector currents arise from the isospin- 
independent and isospin-dependent central and momentum-dependent interactions. These 
currents are short-ranged and numerically far less important than the 7r-like current. 

The model-dependent currents are purely transverse and therefore cannot be directly 
linked to the underlying two-nucleon interaction. The present calculation includes the 
isoscalar pify and isovector unry transition currents as well as the isovector current asso- 
ciated with excitation of intermediate A-isobar resonances (for the values of the various 
coupling constants and cutoff masses in the monopole form factors at the meson-baryon 
vertices, see Ref. pl|). Among the model-dependent currents, those associated with the 
A-isobar are the most important ones. In the present calculation, these currents are 
treated within the static A approximation. While this is sufficiently accurate for our 
purposes here, it is important to realize that such an approach can lead to a gross overes- 
timate of A contributions in electro-weak transitions (see Refs. f22] , f23[ , 24] for a discussion 



of this issue within the context of neutron and proton radiative captures on deuteron and 
3 He, and the proton weak capture on 3 He). 

Finally, it is worth pointing out that the contributions associated with the pn^y, uTc-f 
and A-excitation mechanisms are, in the regime of low to moderate momentum-transfer 
values of interest here (q < 2 fm _1 ), typically much smaller than those due to the leading 
model-independent 7r-like current [2f 



3.1.2 Two-body charge operators 

While the main parts of the two-body currents are linked to the form of the two-nucleon 
interaction through the continuity equation, the most important two-body charge opera- 
tors are model-dependent, and should be considered as relativistic corrections. Indeed, a 
consistent calculation of two-body charge effects in nuclei would require the inclusion of 
relativistic effects in both the interaction models and nuclear wave functions. There are 
nevertheless rather clear indications for the relevance of two-body charge operators from 
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the failure of the impulse approximation in predicting the deuteron tensor polarization ob- 
servable HU , and charge form factors of the three- and four-nucleon systems |^5], |27j . The 



model commonly used []T9| includes the 7T-, p-, and cu-meson exchange charge operators 
with both isoscalar and isovector components, as well as the (isoscalar) prey and (isovec- 
tor) charge transition couplings, in addition to the single-nucleon Darwin-Foldy and 
spin-orbit relativistic corrections. The it- and p-meson exchange charge operators are 
constructed from the isospin-dependent spin-spin and tensor interactions (those of the 
Argonne fig here), using the same prescription adopted for the corresponding current 
operators. 

It should be emphasized, however, that for q < 2 fm -1 the contributions due to these 
two-body charge operators are very small when compared to those from the one-body 
operator. 

3.2 Weak Operators 

In the Standard Model the vector part of the neutral weak current is related to the isoscalar 
(S) and isovector (V) components of the electromagnetic current, denoted respectively as 

•7,(T l -7,0" 

Js and Jv > Vla 

j°< CT = -2 sm 2 6 w3 r + (1 - 2 sin 2 M jV , (23) 

and therefore the associated one- and two-body weak charge and current operators are 
easily obtained from those given in the preceding section. 

The axial charge and current operators too have one- and two-body terms. Only the 
axial current 

j 5 (q) = Ei 5,1 (q) + EJS 2 (q) (24) 

i i<j 

enters in the calculation of the asymmetry. The axial charge operator is not needed in 
the present work. The one-body axial current is given, to lowest order in 1/m, by 

J-' 1 (q) = -G'A(|Q 2 |)^e i ^ , (25) 
where the nucleon axial form factor is parametrized as 

G * (m = a + i£W • (26) 

Here g& is the nucleon axial coupling constant, g& = 1.2654, and the cutoff mass Ka is 
taken to be 1 GeV/c 2 , as obtained from an analysis of pion electroproduction data [ 2E \ 
and measurements of the reaction p — > p + n p9[ 



There are relativistic corrections to j 5 ' 1 as well as two-body contributions arising from 
7T-, p-, /wr-exchange mechanisms and A excitation [p4"| |. All these effects, however, are 



neglected in the present study. The reasons for doing so are twofold: firstly, axial current 
contributions to the asymmetry are small, since they are proportional to the electron 
neutral weak coupling gy ~ —0.074 (see Eq. ([?])); secondly, axial contributions from two- 
body operators are expected to be at the ~ 1 % level of those due to the one-body operator 
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in Eq. (p5|). For example, in the proton weak capture on proton at keV energies [|30|]-this 
process is induced by the charge-changing axial weak current-the n, p, pre, and A two- 
body operators increase the predicted one-body cross section by 1.5%. Such an estimate 
is expected to hold up also in the quasi-elastic regime being considered here. 



4 Connection with Chiral Perturbation Theory 

The purpose of this section is to relate the operators presented in Sections fO and V>?2 



to ChPT. We consider only up and down quarks, in which case the QCD Lagrangian 
has an approximate SU(2)l x SU(2)r chiral symmetry. This symmetry is spontaneously 
broken by the vacuum to the diagonal SU(2) V subgroup, and three pseudoscalar Gold- 
stone bosons, the pions n a , appear in the spectrum. Chiral symmetry provides important 
constraints on the description of low-momentum processes involving pions. In particular, 
it allows one to estimate the relative size of various contributions. 

To accomplish this, the most general effective Lagrangian with broken SU(2)l x 
SU(2)r is constructed. This effective Lagrangian includes terms with arbitrary num- 
ber of derivatives and powers of the quark masses, however higher- dimension operators 
are suppressed by inverse powers of the characteristic mass scale of QCD, Mqcd ~ 1 
GeV. Thus pion interactions are determined as a power series in (g, m^) / 'Mqcd, where 
q is the typical external three-momentum. At low energies this is a small number and 
hence only the lowest-order terms are considered in this paper. 

We start with the tttt and nN Lagrangians. Details can be found, for example, in Refs. 
IB, g. We lump the pions in a field U = u 2 = 1 + ir ■ n/f v - n 2 /2f 2 + where f v = 93 
MeV is the pion decay constant. We denote the nucleon field of velocity and spin 
by N. In order to obtain the currents required for the computation of the asymmetry 
we must use the proper covariant derivatives. The covariant derivatives on the pion and 
nucleon field are constructed in terms of external vector and axial- vector fields V M and A^ 
in the usual manner. They are given by 

DJJ = d^U - i(V M + Ay)U + i U{V» - A„), (27) 



D^N = d^N + ^[u\d^u]N - l - 



u\V®+Af)u + u{Vf-AfW 



iV-3iV, (0) iV, (28) 



where the superscripts (0) and (3) denote isoscalar and isovector components. It is con- 
venient to construct also other quantities that transform covariantly. For example, 

CV = i[w t ,£> M u]+, (29) 
fl = u^F*u + uF^u\ (30) 

where F^ L = d^ L - d u F^ L - i[F^ L , F^ L ] with F* = V„ + A fl and F£ = V M - A^. 

One can find the relation between the external fields and Z° or photon by considering 
the covariant derivative on the quark fields (see, e.g., Ref. f3"T| ): 



A »=2^ T i z >" (31) 
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(J 



,T. 



— - 2 sin 9 w Q g )Z IM + eQ q A 



i-i- 



2 cos 8w 2 

where Q q is the quark-charge matrix. 

From these building blocks we can write the chiral Lagrangians 

C nn = £Tr[DrrfD»U) + . . . , 



(32) 



(33) 



C nN = rf[iv D + g A S- a}N + —rf 

2m 



[v ■ Df - D ■ D - i g A [S ■ D, vv] 



N 



■— N^[S»,S U ] 
Am 



1 



N+... 



(34) 



where k s = —0.12 and k v = 3.71, and . . . denote terms with more derivatives and/or 
powers of the pion mass. One can also write down interactions containing four or more 



nucleon fields || [S^], which are important for a fully consistent description of systems 
involving two or more nucleons. One and two-body currents can be obtained from these 
interactions. 



4.1 Ordering 

The symmetries allow an infinite number of interactions, so an ordering scheme is nec- 
essary for predictive power. We want to estimate the size of matrix elements of one- or 
two-body currents between NN wave functions. These matrix elements involve: the final 
iViV wave function, a two-nucleon propagator, the current operator, another two-nucleon 
propagator, and the initial iViV wave function. Each two-nucleon propagator brings a 
factor m/q 2 , and each loop a g 3 /(47r) 2 . Apart from loops in the operators themselves, 
matrix elements of one- and two-body operators involve, respectively, one and two loops. 

Let us first consider only strong interactions between nucleons pi. These start at 
O(q /f%), and most structures that appear up to 0(q 3 / / 2 Mq CD ) |32| are present in 



the Argonne fig potential JT2|. Use of phenomenological potentials in conjunction with 



currents derived in ChPT has already proven to be a successful approach 0. 

Contributions to the amplitude M 7 start at 0(e 2 m/q 3 ) with the tree-level one-body 
charge operator. First corrections come in tree-level one-body currents from 0(q/m) 
magnetic corrections in the Lagrangian. Second corrections are of two types: (i) one- 
loop corrections and 0(q 2 /Mq CD ) interactions in one-body currents; and (ii) tree-level 
two-body currents. And so on. Contributions to the amplitude M z follow the same 
pattern, but have an extra overall factor of (9(g 2 /M|) = OiG^q 2 /e 2 ). The leading one 
and two-body contributions are described in Sec. [12] below. 

Electromagnetic interactions between nucleons are clearly higher-order effects, also 
included in the t> 18 potential. Interesting are also the contributions to the asymmetry 
that come not from the exchange of a Z° between electron and deuteron, but from Z° 
exchange between hadrons (and photon exchange between electron and deuteron). Let us 
denote these contributions by M l( - Z \ How does M 7( ^ compare to M z ? 
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Figure 2: Electromagnetic one-body diagrams. Solid lines denote nucleons, and the blob 
represents the one-body current. 



Direct Z exchange between nucleons is an 0{Gp) effect, so suppressed by O^G^f 2 ) 
compared to the strong interaction. On the other hand, the electron-deuteron interaction 
is (D(e 2 /q 2 G fl ) larger in M 7 ^ than in M z . As a consequence, the contribution from Z 
exchange between nucleons to M 7 ^ is 0(e 2 f 2 /q 2 ) compared to the leading contribution 
to M z . Since we are interested in q ~ f n , this is ~ e 2 , a small effect. Yet, we are here 
searching for subleading contributions. As we have seen, two-body currents in M z are 
suppressed by 0(q 2 / Mq CD ) compared to the leading one-body effect. Thus, compared 
to two-body currents in M z , M 7 ^ can be 0(e 2 f 2 MQ CD / 'g 4 ); for momenta q ~ f n , 
this is ~ 1. There exist also potentially larger contributions to M^ z \ The size of the 
parity- violating pion-nucleon coupling might be as large as (9(G m / 2 Mqcd/ '<?) compared 
to the parity-conserving coupling ||33|| . In this case pion exchange in the one-body current 
(anapole form factor), in the potential and in the two-body current could all produce an 
effect 0(Mqcd/q) compared to two-nucleon effects in direct Z° exchange. 

Our calculation involves all contributions suppressed by 0(q 2 /Mq CD ) (and also some 
even more suppressed) compared to leading, except for M' 7 ^. However, contributions 
to M 7 ^ from the parity- violating potential were already examined in Ref. 0. Also, 
although not included in any deuteron calculation, the anapole form factor of the nucleon 
seems too small [33] to be relevant here, and it is reasonable to expect that two-body 
currents stemming from the parity-violating pion-nucleon coupling will be equally unim- 
portant. (This should be eventually confirmed in an explicit calculation.) On the other 
hand, the leading two-body currents in the amplitude M z stemming from pion exchange 
are being calculated here for the first time. 

4.2 One- and two-body currents 

The one-body currents contributing to the processes in Fig. [l|, depicted in Figs. || and [|, 
are given to 0(eq/M QC D) by 



2 



£(<?! + gWiW* - r *) (35) 
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Figure 3: Weak neutral one-body diagrams. Symbols as in Fig. [| 



— v 

Am V 



P 5 (r) 



j 5 (r) 



i 



(36) 

(37) 
(38) 



where a = 7 or 0, the coupling constants g% an d fff'^ are given in Table 0, and Qa is 
the nucleon axial coupling constant. These results are in agreement with those presented 
in Ref. ||]. These coupling constants acquire, in higher orders, a Q 2 dependence, see Ref. 
3~4| . In our calculation we use a phenomenological parametrization of the Q 2 dependence, 
as described in Sees. IO and 1373. 



Form Factor 


7 


Z 


9e 


1 


—2 sin 2 8w 


9 V e 


1 


1 - 2 sin 2 9 W 


9m 


1 + K s 


— (1 + k s ) sin 2 9 W 


9m 


1 + K v 


(l + ^)(l-2sin 2 6 W ) 



Table 1: Coupling constants appearing in one-body currents to 0(eq/ Mq C d) in ChPT. 

Note that C nn contains a n-Z mixing term in the form Z ■ dir°, which turns out to 
be proportional to Q^, the four-momentum transfer. Its contraction with the leptonic 
current produces a contribution proportional to the mass of the electron. In the ex- 
treme relativistic limit for the electron under consideration here, this contribution can be 
neglected. 

The two-body contributions to the processes shown in Fig. [I] are depicted in Figs. |] 
and||, where again contributions from tt-Z° mixing are neglected in the extreme relativistic 
limit. To 0(eq 2 /Mq CD ), they are given in momentum space by 
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Figure 4: Electromagnetic two-body diagrams to 0(eq 2 /Mq CD ) in ChPT. Solid (dashed) 
lines denote nucleons (pions). 



Figure 5: Weak neutral two-body diagrams to 0(eq 2 /Mq CD ) in ChPT. Symbols as in Fig. 



j a (ki,k 2 



3i</s(Ti X T 2 ) z 



v n (k 2 )er 2 ■ k 2 cr 1 - v 7T (h)er 1 ■ k 1 cr 2 
v^(k 2 ) - v w (ki) 



k\ - k 2 2 



(ki - k 2 )cr! ■ ki<r 2 ■ k 2 



(39) 



where kj = — pj with pj (p^) denoting the initial (final) momentum of nucleon i. In 
the formula above, 

3U/J mj + k 2 ' 



(40) 



Note that there is no contribution to the axial current operator or the electromagnetic 
charge operator. There is a contribution to the axial charge operator which, however, 
does not enter in the asymmetry computation. These results are in agreement with Ref. 
fl8l , where the Fourier transform of the above expressions are also given in detail. 

In higher order other currents appear. There exist shorter-range currents, which are 
expected to be smaller than the ones from pion exchange with leading order interactions. 
These higher-order effects are parametrized in our calculation through the Riska prescrip- 
tion as outlined in Sections [D] and |3~2} In particular, we do not use Eq. ( f4"0"D for v n (k), 
but the pseudoscalar component of the vi 8 potential. 
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5 Calculation 



In this section we describe the calculation of the deuteron response functions given in 
Eqs. (fn)-(inj). The deuteron wave function is written as 



\d,M d ) 



u ( r ) v M, , w{r) Md 

„ ^011 „ ^211 



(41) 



where the y L gj are standard spin-angle functions, Xm t is a two-nucleon T, My isospin 
state, and u(r) and iw(r) are the S- and D-wave radial functions. 

In the 2 H(e, e')pn reaction the final state is in the continuum, and its wave function 
is written as 

|q;p,5M 5 ,TM T ) =e^ R ^ MsiTA/T (r) , (42) 

where r = ri — r2 and R = (ri + r2)/2 are the relative and center-of-mass coordinates. 
The incoming-wave scattering-state wave function of the two nucleons having relative 
momentum p and spin-isospin states SM$,TMt is approximated as 



^p,SM s ,TAI T ( r ) 



1 



3 ipr 



-l) 5+T e- ip - r 



1 ,(-) 



g rp All ^—-^ 

Xm s Xm t + —f= Z^ 1 °LST 

V 1 JMj LL> 

ax 



u { L ,[(r;p,JST) - 5 L >L jdpr) 



y^ J SjX T M T , (43) 



where 



olst 

ryJMj / * \ 

^LSMs \P) 



1 - 

E 

M L 



-1 



{LMl, SMs\JMj)Ylml{p) 



(44) 
(45) 



The Slst factor ensures the antisymmetry of the wave function, while the Clebsch-Gordan 
coefficients restrict the sum over L and V . The radial functions u\j L are obtained by 
solving the Schrodinger equation in the J ST channel, and behave asymptotically as 



-u { L ,l(r; p, JST) 



1 



SvLh { t\pr) + {Si IL 



(46) 



where Sjf^ is the S"-matrix in the JST channel and the Hankel functions are defined as 
h>L \ x ) = 3l{ x ) ± ji and n L being the spherical Bessel and Neumann functions, 

respectively. In the absence of interactions, u^' L (r;p, JST)/r — > ^i'iii(p r )> an d ip^~\r) 
reduces to an antisymmetric plane wave. Interactions effects are retained in all partial 
waves with J < J max - In the quasi-elastic regime of interest here, it is found that these 
interaction effects are negligible for J max > 7. 

The response functions are written as (only i?^' 7 is given below for illustration) 
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RY r (?,w)= E i# 7 (?,w;£,T), (47) 

S,T=0,1 



where the contributions from the individual spin-isospin states are 
Rr(q,u;S,T) = ~ £ [j^hAl^vMsM^L + E,- 

3 M d M s J ( 27r ) 2 V 




6m m 



with Ag T defined as 



Al T (q, p; M s M d ) = (q; p, SM S T, M T = 0|^(q)|d, M d > . (49) 

Here ^ = —2.225 MeV is the deuteron ground-state energy, the factor 1/2 in Eq. (|H8D is 
included to avoid double counting, and the states |d, MA and |q; p, SMsT, Mt = 0) are 
represented by the wave functions in Eqs. ( ]4"ID and fl4"3"|), respectively. By integrating out 
the energy-conserving 5-function one finds: 

i^( g ,u,;S,T) = -^ 2 r 1 rf(cos^ p )|^ T (g,p,cos^;M 5 M (i )| 2 , (50) 

48 ^ M d ,Ms 1 



where the magnitude of the relative momentum p is fixed by p = y ra(u; + i^) — g 2 / 4, and 
6 p is the angle between q and p. The initial- and final-state wave functions are written as 
vectors in the spin-isospin space of the two nucleons for any given spatial configuration r. 
For the given r the state vector p 7 (q)|<i, Ma) is calculated with the same methods used in 
quantum Monte Carlo calculations of, for example, the charge and magnetic form factors 



of the trinucleons p5| . The r and 9 p integrations required to calculate the amplitudes 
and response function are then performed by means of Gaussian quadratures. 
Finally, note that, since the deuteron is a T = state, one finds 

Rl\q,U]S,T = 0) = -2sin 2 6 w RY{q,u; S,T = 0) , (51) 
R2°(q,u;S,T=l) = (l-2sm 2 6 w )Rr(q,^S,T=l) , (52) 

with similar relations holding between the transverse response functions. 



6 Results and Conclusions 

The asymmetry has been calculated at the kinematics relevant to the SAMPLE exper- 
iment. The incident electron energy was set to E = 193 MeV. SAMPLE measures the 
asymmetry at four different angles (8 = 138.4°, 145.9°, 154.0°, 160.4°). Different electron 
final energies E' correspond to different momentum transfers Q 2 , \Q 2 \ ~ 0.1 GeV 2 in 
the SAMPLE experiment, which is small enough to justify the use of a non-relativistic 
formalism with leading interactions obtained from ChPT. 

The calculated asymmetries, as functions of the electron final energy, are shown in Figs. 
H 0, §, and H for the four different electron scattering angles. For each set of kinematics, 
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the left panels display the asymmetry and the total inclusive cross section, with different 
curves representing one-body contributions, one- plus two-body contributions from pion- 
exchange currents only, and the sum of all contributions. The ratios of one- plus two- 
body contributions from pion only and full currents to one-body contributions for both 
asymmetries and cross sections are displayed in the right panels. 

As is apparent from the figures, the results at all angles are qualitatively similar. Near 
the quasi-elastic peak two-body effects in the asymmetry are negligible, less than 1%, while 
away from the quasi-elastic peak they become relatively more important, increasing the 
asymmetry by at the most 3%. Note, however, that the two-body current contributions 
are large in the inclusive cross section, indeed dominant in the left-hand side of the 
quasielastic peak. In this region the contribution associated with the currents of pion 
range is more than 50% of the total two-body contribution. 

It is interesting to examine more closely the reasons for the relative unimportance of 
two-body current contributions in the asymmetry. At backward angles, the expression for 
the asymmetry can be approximated as 

A R 1 / + (-1 + 4sm 2 6 w )R^ 

G„QV(2>/5a) ~ R¥ ' { ' 

where terms proportional to the longitudinal response functions are suppressed by the 
factor vl/vt < l/tan 2 (#/2), a small number at the angles under consideration here (note 
that the i?J a and R^ a response functions are of the same order of magnitude). It is useful 
to identify the contributions from T = and T = 1 pn final states, and to use Eqs. (|5T|) 
and (j5J), relating the TV/{T = 0, 1) to R^(T = 0, 1). One then finds: 



A 1 — 2 sin 2 6 w {l + r^) + (-1 +4sin 2 #vi/)r 7 ' 5 



G / ,g 2 /(2 v / 2a) 1 + 

r 7 ' 7 = i?W(T = 0)/R^(T = 1) , (55) 
r 7 > 5 = R% 5 (T = l)/i? 7 ' 7 (T = 1) . (56) 

Note that the R^ 5 response function only receives contributions from T = 1 pn final 
states, since the current j 5 ' 1 is isovector. The ratio r 7 ' 7 is much smaller than one, since 
the transverse response is predominantly isovector. For example, at E' = 55 MeV and 
9 = 160.5°, R% J (T = 0) = 0.769 x 10" 5 (0.935 x 10" 5 ) MeV" 1 and R^(T = 1) = 
10.3 x 10~ 5 (25.9 x 10~ 5 ) MeV" 1 , and hence r 7 ' 7 = 0.0748 (0.0361) with one-body (full) 
currents. In contrast, the ratio r 7 ' 5 is of order one; again at E' = 55 MeV and 9 = 160.5°, 
i? 7 ' 5 (T = 1) = -18.3 x 10~ 5 (-26.0 x 10~ 5 ) MeV" 1 , and hence r 7 ' 5 = -1.78 (-1.00) with 
one-body (full) currents. However, it is multiplied by the small factor (—1 +4sm 2 9w) = 
—0.074, and so the asymmetry turns out to be largely independent of nuclear structure 
details. 

Finally, if A denotes the asymmetry obtained by ignoring the contribution of the 
axial current, one finds 

\A\ (-1 + Asm 2 e w )r^ 

~ 1 + 1 o„:„2/i /i , 37^ • ( 57 ) 



\A°\ 1 — 2 sin 6w(l + r7. 
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The computed value for this ratio is shown in Fig. |10| for one of the kinematics of the 



SAMPLE experiment. The contribution of the axial current to the asymmetry is of the 



order of 13% to 24% throughout the kinematical range considered. Note that in Fig. 10 
we have included the small contribution from the longitudinal response. 

As a last remark, we should emphasize that the calculated transverse (R^ 1 ) and 
longitudinal (-RJ 7 ) response functions -including one- and two-body operators- reproduce 
existing Bates data |]37] . 



Since we have performed the computations at the SAMPLE kinematics, the above 
results may be used to account for two-body current corrections in the analysis of the 
experimental data. The SAMPLE experiment measures a convolution of the asymmetry 
A and the cross section (da/dQdE') over a certain range of electron final energies: 

_ JAjda/dndE') dE> 
^totai- !{da/dndE , )dE , ■ W 

The goal of the experiment is to extract the one-body part of A tota i ■ To accomplish this, a 
model that includes only one-body contributions is used to generate the cross section. One 
can now use our results for the ratios of total to one-body contributions in the asymmetry 
and cross section to adjust for two-body effects in the experiment. 

In conclusion, we have presented a fairly complete calculation of the asymmetry in 
quasi-elastic electron-deuteron scattering arising from Z° exchange. Since we find that, 
when the cross section is large at the quasi-elastic peak, the change in the asymmetry 
due to two-body currents is negligible, we expect that these two-body corrections will 
produce a modification in the analysis of the SAMPLE experiment at the % level, too 
small to affect significantly the extraction of the strange and axial form factors of the 
nucleon. It remains to be examined whether the same holds for effects from Z° exchange 
that manifest themselves within the two-nucleon system through the parity- violating pion- 
nucleon coupling. Work along these lines is in progress. 
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Figure 6: Results for scattering of an electron with incident energy E = 193 MeV on 
a deuteron at rest, as function of the electron final energy E' in MeV, for a scattering 
angle 9 = 160.5°. Left panels: longitudinal asymmetry \A\ (top) and cross section o in 
fm 2 /MeV/sr (bottom). Shown are one-body contributions (dotted line), one- plus two- 
body contributions from pion-exchange currents only (dashed line), and the sum of all 
contributions (solid line). Right panels: ratios of one- plus two-body contributions from 
pion only (dashed line) and full currents (solid line) to one-body contributions for the 
asymmetry (top) and cross section aja\ (bottom). 
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Figure 7: Same as Fig. |, but for 6 = 154.0°. 
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Figure 8: Same as Fig. |, but for 9 = 145.9°. 
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Figure 9: Same as Fig. [], but for 9 = 138.4°. 
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Figure 10: The ratio where \A\ is the full asymmetry and \A°\ is the asymmetry 

without the axial contribution for scattering of an electron with incident energy E = 193 
MeV on a deuteron at rest, as function of the electron final energy E' in MeV, for a 
scattering angle 9 = 160.5°. 
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